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ABSTRACT 


The effects of shoaling inner shelf internal tidal bores and solitons on the bottom 
boundary layer have been observed and analyzed during October 2012 at Monterey Bay 
in Monterey, CA. This research utilized measurements from two different observing 
systems: the long range propeller-driven AUV called Daphne, and a bottom boundary 
layer frame measuring changes in temperature, conductivity, pressure, optics, and 
turbulence during a four week time period. The data collected from these platforms were 
used to capture the strength of the internal tides and their influence on the bottom 
boundary layer as they propagate and entrain sediment up into the water column. The 
focus of this field study, sponsored by the National Science Foundation, was to quantify 
turbulent intensity levels for active events near the bed and determine correlations with 
turbulence higher in the water column and the generation of internal nepheloid layers. 
Typical s dissipation rate values for high turbulent events were le-7 nr/s , and as a result, 
turbidity levels near the bed increased on average by a third. Observations show how 
strong changes in the water column currents, turbidity, and turbulence potentially have a 
large effect on divers and subsurface operating vehicles. Further research is required to 
provide a predictive capability of this variability. 
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I. INTRODUCTION 


A. COASTAL OCEANS 

Vigorous productivity in the coastal oceans can be li nk ed to areas of strong costal 
upwelling. Rich nutrient coastal waters provide habitats for more than 90% of all marine 
organisms (PISCO 2009). Regions of upwelling along the coast are essential to 
maintaining ocean productivity due to the injection of nutrients into the upper water 
column (Cazenave et al. 2011). In addition to strong upwelling, coastal oceans are 
constantly affected by tides, winds, and stratification. Combinations of these forcing 
mechanisms acting on the water column play a significant role in the generation of 
internal solitary waves (ISWs) within coastal oceans. McPhee-Shaw (2007) study shows 
how shoaling internal waves are also a vital instrument for transporting nutrients to inner 
shelf kelp ecosystems. As Internal Tidal Bores (ITBs) propagate on shore across a 
continental shelf into shallower waters it interacts with the sloping topography. Steep 
changes in topography increase near bed shear and slow down the shoaling wave causing 
the wave itself to steepen resulting in nonlinear effects including the production of 
leading solitons ahead of the internal tidal bore. Nonlinear internal solitary waves 
influence numerous scientific areas containing cross-shore transport along the shelf, 
ocean acoustic propagation, coastal engineering, and biological oceanography (Farmer 
and Anni 1999). 

1. Monterey Bay 

Monterey Bay is positioned along the Central California Coast between 35° - 37° 
N Latitude and 121° - 123° W Longitude (Figure 1). The Bay is a major location for 
examining and studying internal tide activity due to the underwater topography and 
coastal ocean processes. Typical continental shelves are defined as the area of the ocean 
floor between 30 and 200 meters deep. The continental shelf of Monterey Bay is 
relatively broad for the US West Coast, about 20km wide. A significant portion of the 
shelf s sediment type is composed of mud, fine silt, and clay and their percentages 
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increase toward the shelf break (Figure 2). The shelf break is typically found to be about 
one hundred meters off shore. 

A deep submarine canyon bisects Monterey Bay’s continental shelf. As the 
Monterey Bay Canyon splits through the middle of the continental shelf, it generates steep 
ridges and areas of sloping, irregular bathymetry. The largest slopes within the canyon are 
located at the southern edge where depths change on the order of one thousand meters over 
a distance of five thousand meters. Monterey Bay Canyon terminates within one hundred 
meters of the shore at the eastern end of the bay. Several studies have been conducted to 
determine the effects on internal tides within the canyon (Petruncio et al. 1998). The shelf 
is frequented by energetic non-linear internal waves and is typically a coastal upwelling 
environment from spring until summer (Cazenave et al. 2011). 


Bathymetry of Monterey Bay 



-122.1 -122 -121.9 -121.8 


longitude [deg.] 

Figure 1. Bathymetric map of Monterey Bay. The black line represents the shoreline. 

Elevation scale is shown with blue representing the lowest depths displayed 
and dark red representing sea level. 
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Figure 2. Continental shelf sediment composition in percentage of sand off Western 
California Coast. Monterey Bay is located in the bottom right of the figure 

(From Edwards 2002). 


a. Tides 

The tides in Monterey Bay are mixed, with lunar, M 2 , semi-diurnal being 
the dominant tide. Therefore, the lunar tide is the likely source for tidally driven 
nonlinear internal wave generation. Over the Northern California continental shelf, 
during the 1982 upwelling season as part of the Coastal Ocean Dynamics Experiment 
(CODE), comparisons with the internal wave theory showed the observed vertical and 
horizontal structure of the semidiurnal band horizontal current and temperature 
fluctuations to be consistent with a first baroclinic mode internal wave with a horizontal 
wavelength of 20-30km (Rosenfeld 1990). 

b. Winds 

During the spring and summer months in Monterey Bay, strong north¬ 
westerly winds affect coastal upwelling and advection of the water in the bay (Figure 3). 
Over this time period, maximum wind stress is toward the equator off the coast of 
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California leading to strong seasonal upwelling in Monterey Bay (Hickey 1998). 
However, winter storms disrupt normal wind patterns and lower the mixed layer in the 
water column resulting in less upwelling. Opposing southwesterly winds along the coast 
cause a northward surface current shoreward of the California Current called the 
Davidson Current. The opposing Davidson Current impedes the transport of cold fresh 
water from the north. Consequently, local wind patterns and the direction of the surface 
currents have a significant impact on the stratification in shallow water columns. 



Figure 3. Illustration of strong northwest winds as they blow over the Monterey Bay 
coastal waters. As the surface waters move offshore cold, nutrient-rich water 
moves toward surface (From www.mbari.org). 


c. Stratification 

Stratification varies seasonally within the Monterey Bay. Hickey (1998) 

observed that stratification in California coastal areas are controlled by large-scale 

advection and upwelling of water masses. During periods of upwelling, cool water is 
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pumped onto the continental shelf along the ocean floor resulting in near surface 
stratification. In the spring and summer months, when upwelling is the strongest, a 
predominate two layer system is created in the water column. Strong stratification, 
expressed as large buoyancy frequencies, in Monterey Bay results in an increase of 
internal tide activity and, more specifically, internal bore activity. The buoyancy or 
Brunt-Visailla frequency (N 2 ), 




p dz 


( 1 ) 


is the measure of static stability in the ocean, where g is gravity, p the local density. 
Small N 2 values indicate weak stratification while large N 2 values indicate strong 
stratification. Past studies have shown that near-surface shelf stratification and increased 
maximum buoyancy frequencies exist during late summer months in Monterey Bay 
(Crabbe 2007). These conditions arise because of low to moderate amplitude winds, 
diurnal afternoon sea breeze, few stonns, and strong solar heating of the surface layer 
during late summer. Spatially continuous strong stratification across the shelf is and high 
necessary for the shoreward propagation of the internal tide into the coast. 


B. INTERNAL TIDAL WAVES 

1. Generation and Nonlinear Propagation 

Internal tides are tidally generated internal waves generated as barotropic tidal 
currents flow over sloping topography. Internal tides are commonly observed at 
continental slopes or breaks, straits, and sills within the ocean. Typical internal tides that 
propagate over the continental shelf have current speeds between 5-20 cms' 1 and 
isopycnal displacement amplitudes between 20-50 m (Huthnance 1989; Foreman 1995). 

One of the common observed generation mechanisms for internal tides at the 
continental shelf break is called the lee-wave mechanism (Holloway 1987). During lee- 
wave formation, a depression in the pycnocline is formed when ebbing tidal waters cross 
over a bottom feature or shelf break (Figure 4). As a result, the depressed isopycnals 
generates a standing lee wave of depression. When the flow ceases during the slack tide, 
the wave of depression begins to propagate both onshore and offshore. The offshore 
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propagation of the wave travels into deep waters where it dissipates with little change to 
the water column. On the other hand, as the onshore wave starts to propagate it shoals, 
becomes nonlinear, and begins self-advection. Self-advection causes waves with larger 
amplitudes to move faster. Therefore, the onshore propagating wave gets steeper, often 
fonning a bore-like wave called an internal tidal bore. As the ITB continues to propagate 
shoreward, the subsequent flood tide will assists in the propagation process. Under 
strong stratification conditions, the steep leading edge can generate bore-like shockwaves 
that form a series of depressions in the water column ahead of the ITB. These 
depressions are commonly referred to as internal tidal solitons (Noble et al. 2009). The 
fonnation of solitons depends on the existence of both intrinsic dispersion and 
nonlinearity in the medium (Apel et al. 2007). The magnitude and energy of the soliton 
are increased as the internal tidal bore continues to propagate onshore. In summary, as 
ocean tides ebb and flow over the continental shelf, some of the barotropic tidal kinetic 
energy is converted to baroclinic kinetic energy through the disturbance of the vertical 
density structure in the water column, which may cause internal bores of pulse-like 
density displacements commonly referred to as solitons (Ostrovsky and Stepahanynays 
1989; Stanton and Ostrovsky 1998). 
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a) 


EBB FLOW 


b) 
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Figure 4. Internal tidal bore and soliton generation by the lee wave mechanism, (a) An 
ebb flow across a shelf break causes a depression in the pycnocline. (b) A 
steep-edged shoreward-propagating bore forms during a slack tide, (c) 
Leading edge will continue to steepen due to dispersive effects and the 
internal bore is assisted by the flood tide, (d) The steep leading edge of the 
bore degenerates into a soliton by the process of the dispersive properties of 
internal waves. (From Holloway 1987). 


2. Internal Nepheloid Layers 

Internal nepheloid layers (INLs) are layers of elevated suspended particulate 
material within the water column. INLs are observed commonly near continental slopes 
and shelves, although little is know about their formation mechanism and dispersion 
rates. One hypothesis is that the creation of such layers may fonn by strong baroclinic 
tidal currents that create elevated turbulence along a critical slope and can suspend 
sediments into the water column (McPhee-Shaw et al. 2004). Previous to McPhee- 
Shaw’s laboratory experiment that found intrusions most intense at supercritical sloping 


7 

























conditions, INLs were thought to occur with the critical reflection of semidiurnal internal 
tides. Puig’s (2004) observations further supports McPhee-Shaw’s discovery by finding 
nepheloid layers near critical angles on the northwestern Alboran slope, southwest 
Mediterranean Sea, mainly around depths of 200m to 500m. 

In Monterey Bay, nonlinear internal waves transport sediment offshore 
(Carter et al. 2005). Johnston et al. (2011) observed, through their phase-averaged 
measurements of INLs and currents in Monterey Bay, internal wave-driven westward 
transport of sediment. The study concluded the variable INL locations in depth and 
latitude suggest mesoscale variability alters the conditions at the shallower shelf 
generation site or along the propagation path of the internal tide. Along the Monterey 
Bay Shelf, intrusions can extend 10-20km offshore propagating at a maximum rate of one 
km per 3.5 hours (McPhee-Shaw 2002; Cazenave et al. 2011). 

3. Environmental Impact 

Internal tidal waves and solitons could have a major effect on biogeochemical cycles 
and global climate in coastal oceans. Harmful algal bloom species have a benthic resting 
state and may be the cause of toxic blooms when the donnant cells are carried up into the 
water column through a process called benthic-interior exchange (Rines et al. 2002; 
McManus et al. 2007). Benthic exchange could explain a general correspondence between 
increased levels of Chlorophyll-a and shelf width in major up welling areas (Patti et. al. 2008). 
In addition to harmful algal bloom, internal tidal waves may significantly contribute to 
energy dissipation, and, therefore, ocean mixing and global climate (Baines 1982). 

4. Naval Relevance 

The United States Navy, particularly the submarine and Special Operation Forces 

(SOF) communities, has taken special interest in any field of research that has a direct 

effect on mission success in order to maintain theatre dominance. A majority of 

submarine and SOF operations are executed in the Coastal Ocean environment where 

ISWs are ubiquitous (Figure 5). Critical to the operation of submarines and other 

underwater vehicles is the crew’s ability to utilize acoustic propagation. Consequently, 

internal waves and solitons produce strong “scintillation” in acoustic wave propagation 
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and therefore can hamper the operator’s ability to use underwater acoustics (Kantha et al. 
2000). Apel et al. (2007) also discusses the effect of internal waves and solitons on 
sound wave propagation. Results from their research show ISWs influence propagation 
of acoustic signals in water and may fonn specific conditions for ducting sound waves. 

In addition to underwater vehicle operations, countennining operations are also 
effected by ISWs. Mines can physically be moved or altered during ISW events. Navy 
Divers and Explosive Ordinates Disposal team members diving on a coastal mission can 
have their visibility altered by the presence of INLs. Thus, the United States Navy has a 
high interest in understanding and predicting coastal tide-forced internal wave activity 
that could affect the outcome of their mission. 

The intelligence preparation of the environment is critical for operational success, 
whether the need is halfway around the world or in U.S. coastal waterways, ports, and 
rivers (Sandel 2009). It is essential to gather oceanographic data before conducting naval 
operations in order to provide information on currents, internal waves, and bottom 
geophysical parameters, to deem an area safe to carry out a mission (Dept, of Navy 
2004). A strong understanding of coastal oceanographic processes can better equip naval 
officers in the implementation of naval objectives, and therefore can significantly alter 
the outcome of a mission. 



Figure 5. Map of the world showing sites where internal solitons have been reported 

(From Apel 2007). 
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II. FIELD EXPERIMENT 


The Naval Postgraduate (NPS) Ocean Turbulence Group, Moss Landing Marine 
Lab and Monterey Bay Aquarium Research Institute (MBARI) combined their efforts 
and technology in order to collect data over the continental shelf in Monterey Bay in a 
series of three experiments. The observations were collected as part of the National 
Science Foundation (NSF) study on benthic exchange events and near boundary mixing 
along the continental shelf. 

A. EXPERIMENT SITE 

The Monterey Bay shelf is an ideal site for the experimental study of temporal 
and spatial development of internal tidal bores and nepheloid layers. Steep continental 
shelf breaks are one of the main features in the bay that provide a suitable location for 
internal tide generation. Mud and fine sediments make up the bottom type of the bay 
dominantly, which aids in the generation of nepheloid layers. Monterey Bay’s wide shelf 
supports “mud belt” deposits from winter storm runoffs from local rivers and streams. 
These characteristics of the Monterey Bay shelf are similar to those found on shelves 
throughout the world ocean, indicating that knowledge gained from process studies in 
Monterey Bay can be applied widely. 

Observations in this thesis were gathered from two main instrument systems, a 
bottom boundary layer (BBL) frame and a long-range autonomous underwater vehicle 
(LRAUV), Daphne, from a process study on the southern Monterey Bay Shelf during the 
month of October 2012. The two systems are described in further detail in the following 
section. The frame was deployed on 1 OCT 2012 at 36° 42’ 511” N and 121° 51’ 69” W 
in a water depth of 63 meters (Figure 6). On 11 OCT 2012, Daphne was deployed and 
operated along a line between 121° 49’ 47” - 121° 51’ 69” W longitude along the 36° 42’ 
57” N latitudinal line. The frame was recovered after twenty days of data collection on 
21 OCT 2012 while the Daphne remained in the water until 30 OCT 2012. During the 
month long data collection, Daphne and the frame were both operational for a ten-day 
period from 11 OCT-21 OCT 2012. 
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Figure 6. Map of the Southern Monterey Bay field site, including the deployment 
locations of the two instruments used in the experiment. The red circle 
represents the frame’s location and the red line shows Daphne’s course. 


B. INSTRUMENTS 

1. LRAUV Daphne 

Over the past several years, engineers in J. Bellingham’s group at MBARI have 
developed and advanced a new class of long range AUVs to support physical, chemical, 
and biological sensing missions. Such missions have the endurance capability of 
spanning over 1,000 kilometers or more (MBARI 2013). The LRAUV operate at speeds 
around one meter per second and with minimal sensors allowing the vehicle to range over 
several thousand kilometers. The long endurance capabilities of MBARI’s new vehicle, 
Daphne (Figure 7), allows observations over the entire continental shelf for a period of 
ten to twenty days, which covers over half of the fortnightly tidal cycle. Daphne’s dead 
reckoning navigation is effective at maintaining a set transect in combination with 
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transect end-point surface GPS fixes. Daphne corrects any deviations to the intended 
course when the LRAUV surfaces and receives its current GPS fix. 

Daphne is relatively easy to handle, about 120 kilograms dry weight, and 
designed to be simple enough for operations by individual laboratories. The forward bow 
section of Daphne is free-flooded and houses the scientific payload including the 
Conductivity Temperature Depth (CTD) and three optical sensors. The CTD provides 
regional descriptions of the oceanographic area’s temperature, and salinity, while the 
LRAUV transits through the water column. Fluorometers, one of the three optical 
sensors, are sensitive to biological backscatter and therefore detect the presence of 
chlorophyll in the water column. In addition to the fluorometer, there are 470 nm and 
650 nm optical backscatter sensors that measure particle suspension of benthic material 
near the seabed. The data analyzed from the three optic sensors are used to generate 
high-resolution time series of the water column optical properties. Furthennore, the 
sensors are used to classify when and where on the continental shelf a nepheloid layer 
fonnation was present. Daphne provides the ability to map boundary layer detachment in 
cross-shelf sections, resolving scales from bottom boundary layer turbulent eddies to full- 
water column depth, and width of the shelf (McPhee-Shaw 2002). 

For the benthic exchange NSF project, the NPS Turbulence Group attached a 
custom designed and fabricated turbulence package at the front of the LRAUV (Figures 
8-9). Included in the turbulence package are microconductivity, microtemperature, and 
Acoustics Current Meter (ACM) transducers. These fast-response instruments are used 
to resolve a wide spatial range of turbulent fluctuations present in the water column. 
Specifically, the microconductivity instrument is constructed of two needles 2mm apart 
from one another. An AC voltage is sent through the needles to measure the AC current 
in the water, to measure conductivity in the water. The microtemperature instrument 
measurers the temperature of the water column directly. Time response for the 
instrument is .007 seconds and has a cutoff frequency of 22.73 Hz. Eight ACM 
transducers are paired at a distance of 20cm on paired rods. Each pair of transducer 
measures a different component of velocity, u, v, or w. In this thesis, the turbulence 
measurements will be used to calculate near-bed dissipation rate of kinetic energy and 
scalar variance. 
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Figure 7. The Daphne AUV that provides new capabilities ranging between propeller 
driven vehicles, which only have the endurance of about a day, and 
buoyancy-driven vehicles (gliders), which can operate for months 
(MBARI photograph). 
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Figure 8. Turbulence package with covers on all the extremely sensitive 
microtemperature, microconductivity, and ACM transducers 
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Figure 9. Fully assembled turbulence package on Daphne with covers removed from the 

sensors. 

2. Bottom Boundary Layer (BBL) Station 

The NPS Turbulence Group, led by Tim Stanton and William Shaw, have 
assembled an instrumented frame (Figure 10) designed to measure stresses, bed 
entrainment processes, and fine sediment fluxes both vertically and laterally at the bottom 
boundary layer of the seafloor. On deployment, the frame was aligned with its long axis 
alongshore to lessen the effects of frame-wake disturbances generated by the cross-shore 
wave motions. Instrumentation on the frame includes 1.2MHz (0.5m vertical resolution) 
and 300KHz (2m vertical resolution) Acoustic Doppler Current Profilers (ADCP), two 
hybrid turbulence probes that measure at 0.5m and 1.5m above the bed, a Bistatic 
Coherent Doppler Velocity and Sediment Profiler (BCDVSP) reporting at 40Hz with 
lcm vertical resolution, and a 25cm-path transmissometer. The turbulence instrument 
packages are designed to directly measure vertical turbulent fluxes and turbulent 
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diffusivity using both scalar-dissipation and eddy-correlation techniques. It is possible to 
use the dissipation technique where turbulence is weak, while the eddy-correlation 
method measure fluxes during strong mixing events near the bed. In this thesis I will 
focus on the dissipation techniques. 

a. Acoustic Doppler Current Profiler (ADCP) 

The ADCP measures current velocity profiles by transmitting sound 
pulses through the water column along four narrow, upward-directed beams. The sound 
scattered off particles within the water column are received and the beam velocities of the 
scatters relative to the ADCP result in Doppler shifts of the returning sound energy to 
each ADCP transducer. Transducer quality is vital for data worth. These transducers 
must be directional and efficient in order to process receive the very weak backscattered 
energy. Along-beam current velocities are then calculated using the Doppler shift. In 
order to obtain the orientation of the transducers, an internal compass and tilt sensors 
were measured in both ADCPs. 

By examining the ADCP’s echo intensity, the effective range for the 
ADCP can be determined. Echo intensity is a measure of the signal strength of the echo 
returning from the ADCP’s transmitted pulse. When a constant low value is found for 
several bins, the backscatter profile indicates that the ADCP noise floor has been reached. 
Beyond this range, the ADCP cannot accurately calculate Doppler shifts (RD Instruments 
1998). The Echo intensity from the ADCP also provides a measure for the density of 
water particles suspended in the water column, which is used in this thesis to estimate 
suspended sediment levels in the interpretation of the LRAUV data. 
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Figure 10. Instrumentation on the BBL Frame contains the following: (A) 1.2MHz (left) 
and 300KHz (right) fast sampling broadband ADCP (B) Hybrid turbulence 
probes located 0.5 and 1.5 meters above bed (C) BCDVSP that reports at 
40Hz, 3-component velocity at 1-cm intervals over lm range (D) 25cm-path 

transmissometer 
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III. DATA PROCESSING 


A. LRAUV DAPHNE PROCESSING 

Once Daphne was pulled out of the water following the field experiment, the 
MBARI team off-loaded the raw data and information stored by Daphne’s instruments 
and converted them into Matlab data files. 

The raw data were then extracted and sorted into profiles, bins, sections, and 
segments in order to generate temporal and spatial views of the water column across the 
shelf. During the field experiment, Daphne traveled between two prescribed longitudinal 
coordinates transiting either nominally offshore or onshore. The onshore transit consisted 
of profiles over the entire water column, as the LRAUV would start at the surface and 
dive down to the bed before returning back to the surface. On the offshore transit, the 
LRAUV dove down to a depth of six meters height above bed (HAB) where it remained 
for 500m offshore distances. At the end of each 500m stretch, Daphne would conduct a 
mini profile up into the water column of approximately twenty meters HAB, to sample 
the near-bed gradients of scalar properties. 

Over the course of each transit the LRAUV was diving, surfacing, or changing 
directions. In order to analyze the water column properties and locate internal tide 
events, it is necessary to break down and classify the LRAUV’s movements, called 
segments, in order to form profiles or timeseries elements for each transect analyze the 
water column structure and locate internal tide events. Large on-shore profiles, called 
sections or transects, conducted by Daphne resolve water column structure. After the 
segmenting process is complete, the data go through a binning program so the depth of 
the collected data in water column can be distinguished. On the other hand, the bottom 
tracking movement during the offshore sections are used to calculate the turbulence 
levels near the bed. 

A low-pass-filtered, temporal derivative of pressure, which was gathered 
throughout the experiment, was used to determine Daphne’s mode of travel and break the 
raw data into segments. The low pass filter was used in order to avoid having small 
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fluctuations in seawater pressure produces erroneous segment breaks (Figure 11). For 
each of the data record segments Daphne was assigned an indexing number or state 
variable to keep track of its movement (Table 1). 

As soon as the data segments were assigned state variables, they could be put 
through the segmenting process. A segment is the portion of data collected while the 
LRAUV transits during a specific state. For example, the first segment is when Daphne 
began its initial decent through the water column, state 1 from table, until it changes 
direction to begin its rise through the water column. For each segment variables such as 
start and stop times, latitude and longitude, and HAB calculations are stored. After the 
program is finished running in Matlab, the output is a mat file containing the indexing 
data for specific segmented states. 

Once the segments program was completed, another script was written in order to 
vertically bin the timeseries data for each segment. Binning the data is done in one-meter 
increments starting from the deepest depth of the LRAUV’s transit to the surface of the 
water column. Profiles can then be constructed by compiling the vertical bins together in 
order to cover the entire water column. After individual profiles are created, they are 
then combined into sections, which represent a complete transect either offshore or 
onshore of the LRAUV. During the deployment, Daphne was interrupted twice for 
recharging and sensor refurbishment, which is seen as the gaps to Figure 12. In all, 
Daphne’s data set consists of eighty-eight sections over the twenty-day deployment. 
These profile transect timeseries define the evolution of the water column and help draw 
focus to specific days of interest with intensified activity. 
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State Variable 


Filter Level of Onshore Transit 



State of the LRAUV during Transit 



Figure 11. Variations in sea water pressure recognized by Daphnes location relative to 
the surface of the water. (Top panel) A low-pass filter pressure used to 
smooth out raw data in order to properly categorize LRAUV movement. 
(Bottom panel) Illustrates the resultant LRAUV’s motion in the water 
column. During state +1 Daphne is diving, at -1 Daphne is ascending, and at 
0 Daphne is changing directions or stationary. 
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AUV Number 
Assignment 

Description 

Primary Direction of 
Travel 

0 

Change of direction or 
unnessary movement 

Either 

1 

Diving in the water 
column over 20m 

On-Shore 

•1 

Surfacing in the water 
column over 20m 

On-Shore 

2 

Near bed transit 

Off-Shore 

3 

Diving in the water 
column less than 20m 

Off-Shore 

•3 

Surfacing in the water 
column less than 20m 

Off-Shore 


Table 1. Segment definition for the LRAUV’s transit through the water column. The 
left most column is the different ‘state’ numbers for the LRAUV. The 
center column is a description of each state and the far right column is the 
direction of Daphne transiting through the water column. 
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Figure 12. Deployment transit of the LRAUV from 21 OCT to 30 OCT 2013. The two 
large gaps in the data are due to maintenance being done on the LRAUV out 

in the bay. 


1. Turbulence Calculations 

After the segmenting and binning of the raw LRAUV data was completed, 
turbulence calculations for s (epsilon), % (chi), and K T (turbulent thermal diffusivity) were 
made for all number two segments. These bottom-tracking segments are assumed to 
provide the best data collection for determining the turbulence near the bed. Ultimately, 
the calculated turbulent values should reveal bottom mixing associated with features such 
as passing internal tidal bores and or solitons. Unfortunately, floating matter in the water 
column such as jellyfish or seaweed can contaminate the high-resolution instruments 
making the data unusable. 
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a. 


Epsilon 


Epsilon represents the turbulent kinetic energy dissipation rate, with units 
of m 2 s' 3 . The definition for s was formulated by interpreting the friction term in the 
momentum equation after multiplying through by its velocity component resulting in a 
term equal to: 


3v<V 

dz 


( 2 ) 


where v is the viscosity and u ’ is the fluctuation velocity along a single component. 
Epsilon represents the loss of turbulent kinetic energy due to friction at small scales. Its 
numerical value is large in the presence of turbulence and can be calculated from the 
inertial subrange part of the wave number spectrum of the ACM vertical velocity 
measurement on the LRAUV. 


The inertial sub range transfer scale is seen when the wavenumber has a 
slope of -5/3. Two estimates were calculated for wave number ranges between l-10cpm 
and 10-80cpm (Figure 13). The equation for the vertical (cross-stream) spectrum, S w (k), 
from McPhee 2008, is manipulated in order to solve for s. 

2 

S,«) = ^a^ r (3) 

k 3 

The 4/3 constant comes from theoretical assumption of homogeneity in isotropic 
turbulence (Batchelor 1967). Alpha, a, is the Kolmogorov constant for the along-stream 
spectrum and has a value of .51 when k is the angular wave number in cycles per meter. 


b. Chi 


Analogously to s, % is defined as the rate of destruction of thermal 
variance in the water. The definition for % was formulated by interpreting the molecular 
diffusion term in the temperature equation after multiplying through by another factor of 
T resulting in a tenn equal to: 


x = k t 


Vr V 
v az y 


( 4 ) 
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where Kt is the turbulent thennal diffusivity of seawater. From the diffusive tenn, % can 
be defined as the rate at which temperature fluctuations are smeared away by molecular 
diffusion. 

For this thesis, the microconductivity instruments determined chi’s 
numerical value by directly resolving the dissipation range of temperature. Dillon and 
Caldwell 1980 describe % as the temperature variance dissipation rate satisfying: 

j = 6x£)xQ 2 x S CCp (5) 

where the constant 6 is due to the assumption of isotropy. D is the molecular diffusivity 
of seawater, 1.4e-7. Omega squared is the multiplication factor for the pre-emphasis 
frequency used to electronically pre-emphasize the microconductivity timeseries to form 
the gradient spectrum. And, the final term in the equation, S cc , is the variance of the 
thermal microstructure gradient spectrum that resolves beyond the dissipation scale 
(Figure 14). 


c. Turbulent Thermal Diffusivity 


The turbulent thermal diffusivity, K T , is effective diffusivity used to 
represent turbulent transport in circulation models that cannot resolve turbulent mixing 
eddies. It can simply be defined by manipulating equation (4) to say turbulent thermal 
diffusivity is equal to x over the local background temperature gradient squared: 


Kt (clT 


X 


clZ> 


( 6 ) 


The local vertical temperature gradient was calculated for each bottom-tracking segment 
over a 6m vertical scale. A Matlab script was written to first locate each of the bottom 
tracking segments and then the average was taken of the two profiles before and after 
each segment (Figure 15). Once the averages were obtained, the four gradients were then 
averaged together in order to estimate the temperature gradient for each bottom segment. 
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Figure 13. Amplification of the inertial sub range transfer scales. The line from the top 
left to the bottom right represents a -5/3 slope. Large-scale values are 
represented by epsilon 1 and small-scale values are represented by epsilon 2. 
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Figure 14. Thermal microstructure gradient spectrum calculated by the 

microconductivity instruments on board Daphne. The variance of this curve 
resolves beyond the dissipation rate for chi calculations. 
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dT / dZ Segment 112 



Figure 15. The local vertical temperature gradients of the four segments surround the 
bottom-tracking segment. Segment 112 is the bottom-tacking segment and is 
given a value of the mean to the four gradients averaged together. 


B. BBL FRAME PROCESSING 

1. ADCP 

On the 21 st of OCT, the frame was recovered from the experimental site. The 
data from the two ADCP instruments on the frame was extracted and processed in 
Matlab. These high-resolution instruments calculated velocity observations, specifically 
the U and V components, at the 60m-depth med-shelf location. After separating the two 
components of velocity, current speed was computed by summing the square of each 
velocity component and then taking the square root. Vertical shear was calculated by 
taking the vector vertical differential of the current. Acoustic backscatter was computed 
using the sonar equation: 
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EI = SL + SV + Con-20\og(R)-2aR (5) 

where El is the echo intensity in (dB). SL is the source level or transmitted power (dB). 
SV is the water-mass volume backscattering strength (dB). Alpha, a, is the absorption 
coefficient (dB/meter) and has a calculated value of .51 for the 300KHz ADCP and 2.1 
for the 1200KEIz ADCP. And R is the distance from the transducer to the depth cell (m). 
The constant, Con, is included because the measurement is relative rather than absolute. 
Each step along the process can be shown in a time series with a color intensity plot in 
Matlab to show both the current and acoustic backscatter profile timeseries for visual 
analysis of the water movement structure (Figure 16-17). 
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Figure 16. (Top panel) The cross-shelf or U component of velocity of the BBL frame. 

East is represented by the positive values in red and West is negative and in 
blue. (Middle panel) The along-shore (N or S) or V component of velocity. 
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Shear at 6m HAB 


North is positive and red and south is negative and blue. (Bottom panel) 
Calculated backscatter using the sonar equation. 
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Figure 17. (Top panel) Current speeds in meters per second with red representing max 
speed and blue represents still water. (Middle panel) Vector differential of 
current represented by shear. (Bottom panel) Time series of shear at 6m 

HAB. 
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IV. RESULTS 


A. SEDIMENT SUSPENSION 

With an extensive amount of data collected during the NSF project in OCT 2012 
by the BBL frame and the LRAUV, the first objective was to look at temporal and spatial 
views of the water column for the concentration of suspended, near-bed material and its 
relation to water column velocities. The ten-day span, YD 285 - 295, when the Daphne 
and Frame were both collecting data was the initial focus for analysis. Throughout the 
data analysis process for this time period it was important to keep track of where the 
LRAUV was in reference to the frame while it transected the continental shelf. Daphne 
was determined to be closest to the frame when it was 5m off the coast. 

There were many different ways to represent turbidity in the BBL. For this 
discussion, the best attempt to confirm the presence of sediment suspension was to analyze 
a time series of a combination of instruments from both platforms. These instruments 
included the temperature and 650mn OBS on the LRAUV along with the frame’s 
1200KHz ADCP. Three of Daphne’s onshore full profile transits spanning most of YD 
293 were chosen. The onshore profiles were analyzed instead of the offshore near-bed 
profiles in order to get a broader view of the entire water column. A temperature time 
series of the water column is visually important to attempt to relate the stratification levels 
observed to the presence of turbidity in the water column. Optical measurements from 
Daphne will generate a spatial view of any benthic material present in the water column. 
The backscatter calculated from the 1200KHz ADCP on the frame is used in comparison 
with the OBS instrument on Daphne to amplify any benthic material near the bed as it 
passes over the frame. The 1200KHz ADCP was used over the 300KHz because the 
higher frequency will resolve benthic material better closer to the bed. 

1. YD 293.04 - 293.12 of Sediment Suspension Event 

While spatially tracking the particulate signature during YD 293.04-293.12 
Daphne captures a great example of an internal nepheloid layer in figure 18, panel 2. In 
this panel, the turbidity seen by Daphne as it travels through the water column is 
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highlighted. A particularly high concentration of BBL turbidity is located around the 
60m-depth range. At this depth, it is clear that the sediment spreads out laterally into the 
interior water column and into deeper water. The sediment remains suspended at its 
original depth level and temperature while propagating off shore, indicating the 
generation of an INL. Another area depicted having high turbidity levels is near the coast 
(right side of the panel) toward the top of the water column. It is speculated that this high 
concentration area is due to high levels of phytoplankton in the water column. 

Panel 1 of figure 18 shows the temperature of the water column as Daphne 
transits onshore. By visual examination, the temperatures across the shelf would suggest 
high stratification levels. There is a well-defined two-layered thermal system present 
throughout the water column. The warmest water is concentrated at the surface of the 
water column and is separated from lower depths by a sharp thermocline. An upwelling 
of relatively cold water onto the mid and inner shelf appears to be necessary to allow the 
internal tidal bores to propagate onshore and create internal nepheloid layers. 

The backscatter plot in panel 3 represents the turbidity seen near the seafloor by 
the 1200KHz ADCP of the frame. The ADCP observes high levels of turbidity 
throughout the first four meters of the water column. Highest turbidity levels are 
observed between YD 293.08-293.09 reaching about 6m HAB. These levels of turbidity 
support the quantities observed by the LRAUV as it passed the frame. Keeping in mind 
where the frame is located in reference to the LRAUV, a black line on panels one and 
two was inserted to represent the positioning of the frame, 5m offshore. 

2. YD 293.44 - 293.53 of Sediment Suspension Event 

A few hours after the first offshore transit, Daphne began its second near bed 
observation of the day. The first panel on figure 19 shows the stratification in the water 
column has remained close to the same from the transit earlier in the day. Cold isotherms 
(10 and 11C) have moved upwards and onshore close to the seabed. No relevant change 
was observed to the warmer water that was close to the surface of the water column. 

Noteworthy changes have occurred in the OBS instrument on Daphne, panel 2, 
since the earlier transit. The INL has spread up the continental shelf into shallower 
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waters and the turbidity has become stronger. As a result, higher levels of turbidity are 
observed near the BBL closer to shore. At a distance of 4km off shore there is a clear 
fonnation of a nepheloid layer. Turbidity levels were found elevated in the water column 
about three orders of magnitude greater than the values of the first transit. Hypothesized 
phytoplankton levels near the surface at the coastline have slightly increased from earlier 
in the day. 

Early reported levels during the time series of the backscatter show little to no 
turbidity present in the BBL, panel 3. As the day progresses, a peak in near-bed turbidity 
passes over the frame’s ADCP where a steep increase is reported in suspended sediment. 
The well-observed turbidity region spread as high as the 7m HAB. The peak in turbidity 
correlates with the large-scale internal motions seen in the LRAUV section at 4km 
offshore. During the remainder of the section, turbidity levels began to decrease but at a 
much slower rate than they appeared. 

3. YD 293.89 - 293.98 of Sediment Suspension Event 

During the final transit of YD 293, panel 1 of Figure 20 shows the warmer waters 
have spread deeper into the water column, resulting in the colder isotherms near the bed 
being pushed offshore. As a result, stratification levels are lower then previously 
observed during the first two transits. Additionally, the distinctive two-layered system 
has become less definite. 

Turbidity, panel 2, over the shelf is highest just off shore in no more than 40m of 
water. Over the rest of the shelf, sediment suspension is still observed but at depths much 
closer to the sea floor. There are no distinct INLs observed in the spatial view of the 
water column during this section. The backscatter plot supports the LRAUV optics 
recordings of small amounts of turbidity found near the seabed. Sediment suspension 
over the frame is only observed in the first 2m HAB. 

4. Velocity and Shear during Sediment Suspension Event 

To provide a context for the LRAUV transect profile and constant elevation 
sampling, the ADCP timeseries from the 60m deep fixed ADCPs has been analyzed. 
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Figure 21 is a time series of the ADCP profile timeseries results spanning from YD 
292.5-294.5. The top two panels in the figure show the u and v components of velocity 
above the 300KHz ADCP. U velocity or cross-shore velocity describes the East 
(red/positive) and West (blue/negative) directions of currents. V velocity or along-shore 
velocity defines the North (red/positive) and South (blue/negative) currents. Both 
components experience directional changes to the water column over the two-day period. 
Around YD 293.5 a strong northerly flow with little cross-shore velocities was observed 
near the bed. The strong northerly flow was during the same time the INL present near 
the frame and indicates its direction of propagation. Near the bed a short time later, the 
along-shore velocity shifts to the south and the on shore velocity becomes dominated by 
westerly flow. 

Current speeds are calculated by taking the absolute value of the u and v velocity 
components. By taking the vertical difference of the vector current profile, panel 4 was 
calculated and displays the shear magnitude in the water column. Shear is generated by 
the force of friction acting against the bed of the continental shelf. In order to better 
illustrate the near-bed shear experienced by the offshore LRAUV transits, a time series of 
the shear at 5m HAB was constructed in panel 5. The results from the shear at 5m HAB 
correlates well with the turbidity levels observed in the backscatter plot in panel 3 with 
the strongest shear observed between YD 293.5 and 293.9, the same time period as when 
the INL was present. 
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Figure 18. Suspended sediment event: YD 293.04-293.12. Solid black vertical line 
represents the positioning of the frame as Daphne transits the shelf. The 
curved black line is the bathymetry of Monterey Bay. (Panel 1) Temperature 
of the water column as Daphne transits onshore. (Panel 2) Turbidity seen by 
Daphne as it propagates through the water column. (Panel 3) turbidity seen 
near the seafloor by the 1200KHz ADCP of the frame. 
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Figure 19. Suspended sediment event: YD 293.44-293.53. Solid black vertical line 
represents the positioning of the frame as Daphne transits the shelf The 
curved black line is the bathymetry of Monterey Bay. (Panel 1) Temperature 
of the water column as Daphne transits onshore. (Panel 2) Turbidity seen by 
Daphne as it propagates through the water column. (Panel 3) turbidity seen 
near the seafloor by the 1200KHz ADCP of the frame. 
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Figure 20. Suspended sediment event: YD 293.89-293.98. Solid black vertical line 
represents the positioning of the frame as Daphne transits the shelf. The 
curved black line is the bathymetry of Monterey Bay. (Panel 1) Temperature 
of the water column as Daphne transits onshore. (Panel 2) Turbidity seen by 
Daphne as it propagates through the water column. (Panel 3) turbidity seen 
near the seafloor by the 1200KHz ADCP of the frame. 
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Shear at 5m HAB 
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Figure 21. BBL Frame instruments covering sediment suspension event YD 293.04- 
293.98. (Panel 1 and 2) Velocity (u and v) magnitudes propagating over 
300KFIz ADCP. (Panel 3) Time series of backscatter from 1200KHz ADCP. 
(Panel 4) Shear values calculated from u and v velocity component of 
300KHz ADCP. (Panel 5) Time series of calculated shear values at 5m 

HAB. 
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B. NEAR BED TURBULENCE 

Previous research, in addition to the results from the previous chapter of this 
thesis, has shown the presence of INLs along the continental shelf of Monterey Bay. The 
answer to how these INL are generated has only been hypothesized in the past due to the 
complexity of these events. Motivation to understand how INL events are generated 
leads this analysis of the turbulence data collected from Daphne. Over the twenty days 
Daphne was collecting data, the turbulence instruments were recording high-resolution 
turbulence data. Only a subset of those twenty days had good microconductivity and 
ACM data for the sections after YD 302. In total, nine near bed sections between YD 
302.1 and 303.77 provided useful turbulence data. Unfortunately, when the turbulence 
instruments were sampling the best data, the frame had already been recovered from the 
water; preventing analyses that could be made in conjunction with ADCP timeseries. In 
this particular analysis, the offshore sections of the LRAUV became the focal point for 
results. Offshore sections were examined by calculating the turbulence parameters for 
the bottom tracking segments. 

In order to best describe internal tidal bores and solitons over the Monterey Bay 
Shelf, the offshore transit data from the LRAUV was broken down into multiple six- 
panel figures. The panels consist of the LRAUV height above bed, temperature, OBS, 
and the three turbulence calculations: s, %, and Kt- In addition, for all bottom-tracking 
segments, the thermal gradient spectra were calculated from the microconductivity 
sensor. For this discussion, ITB events were classified into two distinct subgroups: bore 
events or bore events with solitons. 

1. Internal Tidal Bore Event 

Figure 22 spans YD 303.28-303.36 and is characterized as an internal tidal bore 
event. Such classifications can be made by first analyzing the s values, panel 3, for 
turbulence levels. Between YD 303.29 and 303.32, the near bed segments of section 76, 
highlighted in blue, green, red, and cyan, show exceptionally high turbulence values, 
indicating the passing of an ITB. After observing high turbulence near the bed, it is 
important to then examine the temperature gradients, panel 2. There is not much of a 
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temperature gradient observed during these segments and are fairly consistent, therefore 
confirming that no solitons are present. OBS calculations, panel 3, show quite high sediment 
suspension along the bed throughout the entire transit offshore. Even after the turbulence has 
subsided, indicated by lower s values, the turbidity near the sea floor remains. 

The thennal gradient spectra are presented in Figure 23. As expected, the 
observed segments with high turbulence have much higher thermal gradient spectral 
levels than those with low s values. The spectral peak and rapid roll-off at higher 
frequencies for turbulence is consistent with the Batchelor spectrum, and is evidence of 
high-quality turbulence data. 

2. Soliton Events 

With some of the high-energy ITB events, packets of solitons on the leading edge 
of the bore can be seen. These particular events had a significant impact on the water 
column structure as the LRAUV “flew” through them. An ITB with a soliton on the 
leading edge is an extremely nonlinear baroclinic process. Isotherm displacements and 
increased turbulence levels that change suddenly, rather than sinusoidal oscillations, 
indicate the nonlinear baroclinic process. 

An example of the presence of a soliton is seen during section 64 (Figure 24). By 
examining each near bed segment, the pulse-like changes in temperature at 5m HAB 
indicate the presence of a train of solitons followed by the ITB. During the first segment 
(blue) the water column is fairly quiet, with no isothenn displacement, turbidity, or 
turbulence detected. However, on the second segment (green) a very energetic soliton set 
passes through the BBF, indicated by oscillating isotherms. As a result, turbidity in the 
water column rises and falls with the pulse-like solitons. In addition to the increased 
turbidity, the BBF becomes much more turbulent. After the soliton on the leading edge 
of the ITB has passed, turbulence levels maintain a high value while the rest of the bore 
passes across the shelf. 

Only a day after the first soliton was observed, Daphne measures a second one as 
it was conducting its near bed transits. The first segment of section 74 captures the 
highly energetic soliton. Farge fluctuations in temperature were observed as the FRAUV 
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sampled 5m off the sea floor. Each time an isotherm displaced from warm to colder 
water, an increase in turbidity was observed. Similar to the first soliton event, turbulence 
levels were increased near the bed within the soliton, but the turbulence levels decreased 
after the soliton passed. Therefore, the thermal gradient spectrum only shows the first 
segment energetic with turbulence and the subsequent segment inactive. 

Finally, during one of Daphne’s final offshore transits, a third soliton event 
disturbed the water column (Figure 28). Large oscillations were observed in the 
temperature spatial sampling over the first two segments (blue and green). Depressions 
in the isotherm, increased turbidity, and increased turbulence were all observed during 
the same time along the transit. Following the soliton event, the turbulence levels in the 
water column varied between high and lower levels. The spectrum associated with the 
soliton event produce well-shaped Batchelor spectra that fall off at high frequencies. 
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Figure 22. Offshore transit of Daphne section number 76. (Panel 1) Lateral view of the 
LRAUVs movements throughout the water column with segment numbers 
over the bottom tracking segments. (Panel 2) Temperature observed. (Panel 
3) OBS from Daphne. (Panel 4-6) Turbulence calculations for epsilon, chi, 

and Kj 
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Figure 23. Thermal microstructure gradient spectrum calculated by the 

microconductivity instruments on board Daphne of each bottom-tracking 
segment over section number 76. Colors correlate with previous figure of the 

same offshore transit section. 
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Figure 24. Offshore transit of Daphne section number 64. (Panel 1) Lateral view of the 
LRAUVs movements throughout the water column with segment numbers 
over the bottom tracking segments. (Panel 2) Temperature observed. (Panel 
3) OBS from Daphne. (Panel 4-6) Turbulence calculations for epsilon, chi, 

and Kt 
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Figure 25. Thermal microstructure gradient spectrum calculated by the 

microconductivity instruments on board Daphne of each bottom-tracking 
segment over section number 64. Colors correlate with previous figure of the 

same offshore transit section. 
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Figure 26. Offshore transit of Daphne section number 74. (Panel 1) Lateral view of the 
LRAUVs movements throughout the water column with segment numbers 
over the bottom tracking segments. (Panel 2) Temperature observed. (Panel 
3) OBS from Daphne. (Panel 4-6) Turbulence calculations for epsilon, chi, 

and Kj 
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Figure 27. Thermal microstructure gradient spectrum calculated by the 

microconductivity instruments on board Daphne of each bottom-tracking 
segment over section number 76. Colors correlate with previous figure of the 

same offshore transit. 
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Section 80 
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Figure 28. Offshore transit of Daphne section number 80. (Panel 1) Lateral view of the 
LRAUVs movements throughout the water column with segment numbers 
over the bottom tracking segments. (Panel 2) Temperature observed. (Panel 
3) OBS from Daphne. (Panel 4-6) Turbulence calculations for epsilon, chi, 

and Kt 
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Figure 29. Thermal microstructure gradient spectrum calculated by the 

microconductivity instruments on board Daphne of each bottom-tracking 
segment over section number 80. Colors correlate with previous figure of the 

same offshore transit. 
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V. DISCUSSION 


A. SUSPENDED SEDIMENT EVENTS 

Monterey Bay experiences high levels of near-surface stratification during late 
summer, approximating a two-layered system. Over the course of the NSF experiment, 
high-level stratification events were captured by the temperature profile transects 
measured by the LRAUV as it transited the water column. During the three transects by 
Daphne on YD 293 a transition from a strong two-layered stratification in the morning to 
a less defined and weaker stratification at night was observed. At this time cold 
isotherms were pushed onshore due to the effect of internal tides. Other researchers have 
shown how the processes of coastal upwelling, which raises the stratification levels, 
amplify this shifting in cold isotherms onshore. 

As hypothesized, the origin of the INL was found to be between 60-70m depths 
and subsequently, the suspended sediment was observed to spread horizontally from the 
origin as far as 2.5km along the continental shelf. After suspension has occurred, the INL 
seems to propagate up the shelf. Support for shore ward propagation is found in the high 
concentrated levels of turbidity found near the shore at the end of YD 293. 

By further analyzing the velocity measurements captured from the 300KHz 
ADCP, it can be detennined that during YD 293 an ITB was propagating down from the 
north. During this time, large amounts of shear and turbidity were present near the sea 
floor and influenced the generation of the INL. In order to further explore how the INLs 
are created and propagate, turbulence levels were analyzed. 

B. TURBULENCE EFFECTS 

Internal tidal bores and solitons were frequently observed providing the necessary 
stress and turbulence levels near the bed to suspend sediment into the water column. 
During section 76, high levels of turbulence were observed without any presence of high 
frequency solitons at its leading edge. After the highly turbulent bore passed and 
increased the level of turbidity near the bed, the suspended sediment did not quickly 

return to the seabed. Therefore, turbidity caused by ITBs can remain in the water column 
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for multiple hours after the bore has passed. During active periods of internal tidal waves 
it is difficult to detennine exactly how long the turbidity levels remain suspended because 
currents advected the INLs. 

When solitons were observed around depth levels of the frame, section 64, they 
were closely followed by their accompanied ITB. The results for this section confirm 
that solitons are found on the leading edge of the ITB and propagate in strong 
unidirectional pulses of current over the shelf. During analysis of sections 74 and 80, 
high levels of turbulence did not follow the solitons but rather they were followed by a 
period of low turbulence. One explanation to why the ITB was not seen following the 
soliton can come from the detection location of the soliton on the shelf. After solitons 
have been generated at the leading edge of an ITB, they have the ability to break away 
and propagate at a much faster rate than the ITB. Therefore, when solitons are found 
close to the shoreline and followed by periods of calm turbulence, it can be speculated 
that these highly non-linear soliton pulses have propagated ahead of the ITB. 

ITBs and solitons create enough turbulence within the BBL that high levels of 
turbidity are seen as they pass though. Surface wave forcing can help to suspend and 
entrain the sediment, but these observations suggest that the solitons and tidal bores are 
the dominant force for transport up into the water column and net transport of sediment. 
On the other hand, local high turbulence is not a necessary condition for high levels of 
suspended sediment to exist. For example, in figure 26, the second segment has low 
turbulence values but a rather high level of turbidity. Such events can occur in the water 
column because of advection and the relatively long settling time of the suspended 
sediment. 
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VI. CONCLUSION AND RECOMMENDATION 


A. SUMMARY 

The effects of shoaling inner shelf internal tidal bores and solitons on the bottom 
boundary layer have been observed and analyzed during October 2012 from Moss 
Landing to Monterey Bay in Monterey, CA. This research utilized measurements from 
two different observing systems: the long range propeller-driven AUV called Daphne, 
and a bottom boundary layer frame measuring changes in temperature, conductivity, 
pressure, optics, current profiles, and turbulence during a four week time period. The 
data collected from these platforms were used to capture the strength of the internal tides 
and their influence on the bottom boundary layer as they propagate and entrain sediment 
into the upper water column. Internal tidal bores and soliton velocities had components 
in both cross-shore and alongshore directions. Internal bore events may be categorized 
by whether or not they contained solitons. Due to the nonlinear nature of solitons, no 
consistent pattern, shape, duration, or maximum excursion could be linked between each 
event. However, each soliton train was characterized by their steep isotherm excursions. 
In addition, solitons, when present, are always found propagating in front of the bores, 
generated at the leading edge, but can separate away from the slower moving ITB. 

The main focus for this field study, sponsored by the National Science 
Foundation, was to quantify turbulent stress levels for heightened events near the bed. In 
addition, link any possible correlations with turbulence and the generation of internal 
nepheloid layers. Typical s values for high turbulent events were le-7 m 2 /s 3 , with 
turbidity levels near the bed increased on average by a third. During these events it was 
directly observed that fine bed sediment was transported upward into the entrained water 
column where they can form INLs. Observations show how strong changes in the water 
column currents, turbidity, and turbulence would have a large effect on divers and 
subsurface operating vehicles. Further research is required to provide a predictive 
capability of this variability. 
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B. NAVAL IMPACT 

Results from this research are valuable in better understanding the highly complex 
coastal environment and how bathymetry such as submarine canyons and shelf breaks 
can generate tidal forced waves. Regions of strategic military interest, including the 
South China Sea and straits of Gibraltar, are frequently populated by ITB and soliton 
activity. Energetic internal disturbances like bores and solitons can disrupt and alter 
military operation missions. These missions include but are not limited to mine warfare, 
submarine warfare, and special operation warfare. Increased knowledge of the diverse 
coastal environment would result in improvements to littoral operations, nearshore model 
resolution, and a better understanding of the coastal acoustic environment. 

C. RECOMMENDATIONS FOR FUTURE RESEARCH 

Potential research suggestions from this fieldwork include analyzing the data from 
other instruments on the BBL frame. The BCDVSP is a non-invasive, high-resolution 
acoustic instrument that measures turbulent stress within internal waves and any low- 
frequency current boundary layers. By analyzing such data, turbulent vertical and lateral 
sediment fluxes can be estimated at centimeter resolutions. In addition, to the BCDVSP 
there were temperature sensors ever 2m up to 16m HAB attached to the cable connecting 
the frame to the surface buoy. The temperature sensors would aid in resolving any steep 
temperature gradients generated by internal tidal bores and solitons. Furthermore, the 
strength of ITB events could be further characterized by analyzing the wind and tide 
relationship with the currents over the continental shelf. 

Future deployments would benefit from having both the BBL frame and the 
LRAUV operational in order to collect the necessary data to best capture these complex 
tidal events. Utilizing the array of instruments attached to the two platforms will produce 
enough necessary data to better categorize and explain the extremely nonlinear internal 
events better than ever done before. The implementation of the flux package to the front 
of the LRAUV provided excellent turbulence data and would benefit any future study of 
internal tidal bores and solitons. An increasing understanding of the effects that tidal 
bores have on the generation of INLs will strongly enhance the US Navy’s dominance of 
the undersea domain. Ultimately, this understanding will provide the decision makers of 
the future the necessary knowledge to ensure optimal planning for costal naval missions. 
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